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Abstract 

The coexistence of a homogeneous (Meissner-like) phase of spin-triplet super- 
conductivity and ferromagnetism is investigated within the framework of a phe- 
nomenological model of spin-triplet ferromagnetic superconductors. The results are 
discussed in view of application to metallic ferromagnets as UGe2, ZrZn2, URhGe, 
and Fe. 

1. Introduction 

Experiments at low temperatures and high pressure have indicated a coexistence of fer- 
romagnetism and superconductivity in the metallic compounds UGe2 [11121131, ZrZn 2 jlj, 
URhGe [S] and also in Fe [Hj. In contrast to other superconducting materials (see, e.g., 
Refs. |3 IE]); in these metals the phase transition temperature to the ferromagnetic state 
is higher than the phase transition temperature to the superconducting state and the 
superconductivity not only coexists with the ferromagnetic order but is enhanced by it. 
It is widely accepted [TJ [TU] that this superconductivity can be most naturally under- 
stood as a spin-triplet rather than a spin-singlet pairing phenomenon (see, e.g., Ref. [0]). 
The experiments |Hj on high-pressure crystal modification of Fe, which has a hexago- 
nal closed-packed structure, are also interpreted P in favour of the appearance of same 
unconventional superconductivity. Note, that both vortex and Meissner superconductiv- 
ity phases jE] are found in the high-pressure crystal modification of Fe where the strong 
ferromagnetism of the usual bcc iron crystal probably disappears 10,. 

The phenomenological theory that explains coexistence of ferromagnetism and unconven- 
tional spin-triplet superconductivity of Landau- Ginzburg type was derived [TTJ H2j on the 
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basis of general symmetry group arguments . It describes the possible low-order coupling 
between the superconducting and ferromagnetic order parameters and establishes several 
important features of the superconducting vortex state in the ferromagnetic phase of un- 
conventional ferromagnetic superconductors ^Tj ^] . Both experimental and theoretical 
problems of the ferromagnetism in spin-triplet superconductors seem to be quite different 
from those in conventional (s-wave) ferromagnetic superconductors [3IH]. 

In this letter we investigate the conditions for the occurrence of the homogeneous (Meissner- 
like) phase of coexistence of spin-triplet superconductivity and ferromagnetism. Such a 
phase of coexistence may occur at relatively small magnetization and at zero external mag- 
netic field. Taking in mind this circumstance and using model considered in Refs. ITS] 
we show that the phase transition to the superconducting state in ferromagnetic supercon- 
ductors can be either of first or second order and this depends on the model parameters 
that correspond to the particular substance. Our investigation is based on the mean-field 
approximation as well as on familiar results about the possible phases in nonmagnetic 
superconductors with triplet (p-wave) Cooper pairs jTUHH]. We neglect all anisotropies, 
usually given by the respective additional Landau invariants and gradient terms [9, in the 
Ginzburg-Landau free energy of unconventional superconductors. The reasons is that the 
inclusion of crystal anisotropy is related with lengthy formulae and a multivariant analysis 
which will obscure our main aims and results. Let us emphasize that the present results 
should be valid in the same or modified form when the crystal anisotropy is properly 
taken into account. We have to mention also that there is a formal similarity between the 
phase diagram obtained in our investigation and the phase diagram of certain improper 
ferroelectrics [TH] . 

2. Model 

We consider the Ginzburg-Landau free energy [TT| IT2] F = J d 3 xf(ip,Ai), where 

/ = ^(D*^)(D^) + aM 2 + \\^ + a f M 2 + ^M A + ii,M.^ x if) . (1) 

In Eq. (jTJ), Dj = (V — 2ieAj/hc), and Aj (j = 1,2,3) are the components of the vector 
potential A related with the magnetic induction B = V x A. The complex vector ip = 
(ipi,ip2,i J 3) = {4>j} is the superconducting order parameter, corresponding to the spin- 
triplet Cooper pairing and M. = {M.j} is the magnetization. The coupling constant 
7o = 47r J > is given by the ferromagnetic exchange parameter (J > 0). Coefficients 
a s = a S {T — T s ) and a/ = a/(T — Tf) are expressed by the positive parameters a s and 
acf as well as by the superconducting (T s ) and ferromagnetic (Tf) critical temperatures in 
the decoupled case, when Mipiip ^-interaction is ignored; b > and (3 > as usual. 

We assume that the magnetization Ai is uniform, which is a reliable assumption outside a 
quite close vicinity of the magnetic phase transition but keep the spatial (x—) dependence 
of ip. The reason is that the relevant dependence of tp on x is generated by the diamagnetic 
effects arising from the presence of M. and the external magnetic field H IT2*] rather 
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than from fluctuations of ip (this effect is extremely small and can be safely ignored). 
First term in (JJJ will be still present even for H = because of the diamagnetic effect 
created by the magnetization Ai = B/An > 0. As we shall investigate the conditions for 
the occurrence of the Meissner phase where ip is uniform, the spatial dependence of ip 
and, hence, the first term in r.h.s. of (0) will be neglected. This approximation should 
be valid when the lower critical field H c i(T) is greater than the equilibrium value of the 
magnetization Ai in the phase of coexistence of superconductivity and ferromagnetism. 

One may take advantage of the symmetry of model ([TJ and avoid the consideration 
of equivalent thermodynamic states that occur as a result of the respective continuous 
symmetry breaking at the phase transition point but have no effect on thermodynam- 
ics of the system. We shall assume that the magnetization vector is along the z-axis: 
Ai = (0,0, Ai), where Ai > 0. We find convenient to use the following notations: 
<Pj = b x /%, <pj = ^exp(^), M = P^M, 7 = 7o/(&W /4 , r = a./Vb, t = a f /y/J3, 
<P 2 = (0 2 + </> 2 + 2 ), and 9 = (6 2 -9 1 ). 

We shall not dwell on the metastable and unstable phases described by the model (0) 17J 
as well as on the vortex phase flj corresponding to \B\ > H c \. Rather we consider 
the stable homogeneous phases at zero external magnetic field (H = 0) that are described 
by uniform order parameters Ai and ip. We shall essentially use the condition Tf > T s . 

3. Results and discussion 

The model (JIJ describes three stable homogeneous phases: (1) normal phase ((f) j = M = 0) 
(hereafter referred to as N-phase), (2) ferromagnetic (FM-) phase ((f) j = 0, M 2 = — t > 0), 
and (3) a phase of coexistence of superconductivity and ferromagnetism (FS-phase), given 
by 03 = 0, 9 = 2n(k - 1/4), (k = 0, ±1, ...), 2 = (-r + 7M) > 0, 6M 2 > ( 7 2 - 2t), and 

y=(y-i)M-M 3 . (2) 

It is not difficult to determine the domains of existence and stability of the phases N, FM, 
and FS. Note, that here we use the term "condition of stability" to indicate the necessary 
condition of stability when the respective phase corresponds to a minimum of the free 
energy, i.e., in both cases of stable and metastable states. When a phase corresponds to 
a global minimum of the free energy (a sufficient condition of stability) it will called a 
"stable phase." Thus we easily find the following existence and stability regions: (t > 0, 
r > 0) for the N-phase, (t < 0, r > 7M) - for FM. In order to obtain the same domain 
for FS one should consider Eq. (J2J) together with the additional existence and stability 
conditions corresponding to this phase: 7M > r and 3M 2 > Mq, where M > is defined 
by the r(M) = 0; see Eq. ©. 

Eq. (J2J is shown in Fig. 1 for 7 = 1.2 and t = —0.2. For any 7 > and t, the stable 
FS thermodynamic states are given by r(M) < r m = r(M m ) for M > M m > 0, where 
M m corresponds to the maximum of the function r(M). Functions M m (t) and Mo(t) = 
(-t + 7 2 /2) x / 2 = v / 3Af m (t) are drawn in Fig. 2 for 7 = 1.2. Functions r m (t) = 4M^(t)/7 
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Figure 1: h — 77-/2 as a function of M for 7 = 1.2, and t = 



-0.2. 



for t < 7 2 /2 (the line of circles in Fig. 3) and r e (t) = 7|t| 1//2 for t < (the dotted in 
Fig. 3) define the borderlines of existence and stability of FS. 

In the region on the left of the point B (Fig. 3) with coordinates (— 7 2 /4,7 2 /2), FS satisfies 
the existence condition 7M > r only below the dotted line [r < r e ] . In the domain confined 
between the lines of circles and the dotted line on the left of the point B the stability 
condition for FS is satisfied but the existence condition is broken. The inequality r > r e (t) 
is the stability condition of FM for < (-t) < 7 2 /4. For (-t) > 7 2 /4 the FM phase is 
stable for all r > r e (t). The dotted line on the left of the point B, i.e. for (—t) > 7 2 /4), 
is a line of the second order FM-FS phase transition. On this line the equilibrium order 
parameters are given by (f)j = and M eq = y/\t\. Therefore, the phase diagram for 
(—t) > 7 2 /4 is clarified for any r. When r < the FS phase is stable and is described by 
the function r(M) for M > M , as shown in Fig. 3. 

The part of the t-axis given by r = and t > 7 2 /2 in Fig. 3 is a phase transition line of 
second order that describes the N-FS transition. The same transition for < t < 7 2 /2 is 
represented by the solid line AC which is the equilibrium transition line of a first order 
phase transition. This equilibrium transition curve is given by the function 

r eq (t) = \ [3 7 - ( 7 2 + 16t)) 1/2 ] M eq (t), (3) 

where 

M eq (t) = ±= [ 7 2 - St + 7 (7 2 + 16t) 1/2 ] V2 (4) 
is the equilibrium value (jump) of the magnetization. The order of the N-FS transition 
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Figure 2: M versus t for 7 = 1.2: the dashed line represents M , the solid line represents 
M e q, and the dotted line corresponds to M m . 
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Figure 3: The phase diagram in the plane (t, r) with two tricritical points (A and B) and 
a triple point C; 7 = 1.2. 
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changes at the tricritical point A. 

The domain above the solid line AC and below the line of circles for t > is the region 
of a possible overheating of FS. The domain of overcooling of the N-phase is confined by 
the solid line AC and the axes (t > 0, r > 0). At the triple point C with coordinates [0, 
r eq (0) = 7 2 /4] the phases N, FM, and FS coexist. For t < the straight line 

< q (t) = ^ + \t\, - 7 2 /4<t<0, (5) 

describes the extension of the equilibrium phase transition line of the N-FS first order 
transition to negative values of t. For t < (— 7 2 /4) the equilibrium phase transition FM- 
FS is of second order and is given by the dotted line on the left of the point B (the second 
tricritical point in this phase diagram). Along the first order transition line r* (t), given 
by (jSJ), the equilibrium value of M is M eq = 7/2, which implies an equilibrium order 
parameter jump at the FM-FS transition equal to (7/2 — yjtj). On the dotted line of 
the second order FM-FS transition the equilibrium value of M is equal to that of the FM 
phase (M eq = At the triple point C the phases N, FM, and FS coexist. 

In conclusion, we have demonstrated that the model (1) of ferromagnetic spin-triplet 
superconductors gives a quite complex phase diagram containing three ordered phases, 
two types of phase transitions, and two tricritical points, and a triple point. Further 
considerations of the effect of additional terms in the free energy such as terms of 
the type A^ 2 ^! 2 and/or terms describing the Cooper pair and crystal anisotropy IHj 
may give more information about the shape of the phase diagram outlined in the present 
paper. 
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